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Abstract 

We present the detailed theoretical study of a heterostructure comprising of two coupled ferro- 
magnetic superconducting layers. Our model may be also applicable to the layered superconductors 
with alternating interlayer coupling in a parallel magnetic field. It is demonstrated that such sys- 
tems exhibit a competition between the nonuniform Larkin-Ovchinnikov-Fulde-Ferrel (FFLO) state 
and the tt superconducting state where the sign of the superconducting order parameter is opposite 
in adjacent layers. We determine the complete temperature-field phase diagram. In the case of low 
interlayer coupling we obtain a new tt phase inserted within the FFLO phase and located close to 
the usual tricritical point, whereas for strong interlayer coupling the bilayer in the tt state reveals 
a very high paramagnetic limit and the phenomenon of field-induced superconductivity. 
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I. INTRODUCTION 



The question of coexistence of singlet superconductivity and magnetism has been adressed 
for many years. It was found that the superconducting order parameter is destroyed by a 
magnetic field both via the orbital effect [1] and the paramagnetic effect. [2] In the usual case 
of an isotropic three-dimensional (3D) superconductor under an external magnetic field, 
the orbital effect prevails and leads to the well-known temperature-field phase diagram of 
conventional type I or II superconductors. [3] In contrast, superconductivity is essentially 
suppressed by the paramagnetic effect in presence of a ferromagnetic exchange interaction. 
This is also true for quasi-two-dimensional (2D) superconductors under in-plane magnetic 
field and for heavy fermions materials wherein the orbital effect is partially quenched. In 
the whole paper, the magnetism is characterized by an internal exchange field h (given 
in energy units) which may arise either from an externally applied magnetic field or from 
ferromagnetic ordering. Note that ferromagnetism must be weak in order to avoid complete 
suppression of superconductivity. This is realized in rare-earth metals or actinides in which 
the indirect exchange interaction leads to Curie temperatures of a few degrees. 

Superconductors with internal homogeneous exchange field h exhibit a very special be- 
haviour. According to Chandrasekhar[4] and Clogston,[5] at zero temperature uniform su- 
perconductivity should be destroyed when the polarization energy of the free electron gas 
exceeds the energy gain due to Cooper pairing in the BCS ground state. This criterion 
gives the exchange field h p (T = 0) = Ao/a/2 where the superconductor should undergo 
a first-order transition to the normal state, Ao = 1.76T c o being the zero temperature su- 
perconducting gap. Larkin and Ovchinnikov [6] and Fulde and Ferrell [7] (FFLO) pre- 
dicted the existence of a nonuniform superconducting state with higher critical exchange 
field hf£ LO (T = 0) = 0.755A > h p (T = 0) and second-order transition to the normal state. 
This prediction was made for 3D superconductors. In quasi-2D superconductors the critical 
exchange field of the FFLO state is even higher, namely h^ LO {T = 0) = A ,[8] while in 
quasi-one-dimensional systems there is no paramagnetic limit at all. [9] The appearance of 
the modulated FFLO state is related to the pairing of electrons with opposite spins which 
do not have the opposite momenta anymore due to the Zeeman splitting. From now on 
we focus on the 2D case for which a generic temperature-exchange field phase diagram has 
been established. [8] At low field and temperature, the ground state is characterized by a 
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uniform superconducting order parameter. A tricritical point, located at h* = 1.07T c0 and 
T* = 0.56T c0 , is the meeting point of three transition lines separating the normal metal, 
the uniform and the nonuniform superconductors. At T < 0.56T c0 , the (low-field) uniform 
superconductor is separated from the (high-field) normal metal by a narrow FFLO nonuni- 
form superconducting phase. In contrast, at T > 0.56T c0 , the system undergoes merely a 
second-order phase transition from the uniform superconductor to the normal metal when 
increasing the exchange field. The nonuniform FFLO state is settled in a small region of 
the phase diagram and is very sensitive to impurities, [10, 11] making it difficult to observe 
experimentally. Nevertheless, several evidences of the FFLO state have been obtained re- 
cently in organic superconductors [12, 13] and in heavy fermions compounds, see Martin et 
al. [14] and references therein. 

In the context of organic and high-T c superconductors, layered systems made of conduct- 
ing atomic planes have been extensively studied. [15] In order to investigate the interplay 
of superconductivity and magnetism in such anisotropic systems, [16] Andreev et al. con- 
sidered a periodic array of alternating ferromagnetic and superconducting 2D planes. [17] 
Solving the corresponding Gor'kov equations, these authors established the existence of a 7r 
state wherein each F layer separates superconducting planes with opposite order parameter. 
This is relevant for the ruthenocuprate compound RuSr 2 GdCu20s which comprises Cu02 
superconducting planes and Ru0 2 magnetic planes. [18, 19] A related system is an isolated 
F/S/F trilayer which exhibits the so-called superconducting spin-valve effect. Namely, its 
critical temperature is higher in the antiparallel (AP) orientation of the layers magneti- 
zations than in the parallel (P) orientation both for thick layers[20, 21] and atomic size 
layers. [22, 23] Surprisingly, in the atomic thickness limit, the superconducting gap at zero 
temperature is higher for P orientation of the magnetizations. [22-25] Hence one expects a 
transition from AP to P orientation by cooling the system below a finite crossing temper- 
ature. The recent progress in molecular beam epitaxy[26] enables to fabricate such F/S/F 
trilayer with atomic thicknesses. 

In this paper, we consider a periodic array of SF bilayers. Each bilayer is made of two 
atomic planes coupled by single electron tunneling. Both exchange fields and BCS super- 
conducting pairing are present in each SF plane. The possibility of \ = 71 phase difference 
between the planes inside each bilayer is also taken into account. In the whole paper, we 
assume that the coupling t' between successive bilayers is considerably weaker than the 
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intra-bilayer coupling t. Our study is performed within the framework of the BCS theory 
of s-wave superconductivity. Solving exactly the Gor'kov equations in the limit t'/t — > 0, we 
first derive the critical temperature and the superconducting gap both for parallel (P) and 
antiparallel (AP) orientation of the magnetizations. We show that the critical temperature 
is higher for the AP orientation than for the P orientation whereas it is the opposite for 
the zero temperature gap. We also investigate the interlayer Josephson current in the small 
coupling limit: the current increases as a function of the exchange field for AP orientation 
whereas it is field-independent for the P orientation. Furthermore, we find that for low 
exchange fields and high temperatures, the ground state corresponds to identical supercon- 
ducting order parameters on adjacent layers. For high enough fields and/or low enough 
temperatures, the tt phase ground state is favoured and compete with the FFLO state. For 
the P orientation, the full temperature-exchange field phase diagram is constructed in the 
two limits of extremely low and high coupling between the planes. As expected, for per- 
turbative coupling between two SF planes, the phase diagram is very close to the quasi-2D 
superconductor's phase diagram. Nevertheless an important change arises. Indeed a new 
7i— phase is inserted inside the usual FFLO phase close to the tricritical point. For higher 
tunneling coupling t > T c0 , this n— phase is pushed to low temperatures T < T c 2 /t and high 
fields h ~ t. In this unusual superconducting phase, the Zeeman splitting is compensated 
by the bonding/antibonding energy splitting due to single-electron tunneling between the 
planes. [27] As a result, field-induced superconductivity and enhanced paramagnetic limit 
are realized in this simple model. These new phenomena are encountered due to the intro- 
duction of an additional discrete degree of freedom, here the layer index j. The layer index 
acts as a pseudo-spin and thus enlarges the usual spin-space for singlet pairing. This idea 
was introduced by Kulic and Hofmann[28] in the context of two-bands superconductivity for 
which the pseudo-spin was the band index. Nevertheless, these authors did not investigate 
the presently studied ir state. 

The outline of the paper is the following. In Sec. II, we present the model, derive the 
corresponding Gor'kov equations and give their exact solutions. In Sec. Ill, we investigate 
the critical temperature, the gap and the interlayer Josephson current in the small exchange 
field regime for which there are only uniform superconducting phases. In the last two sections 
the temperature-exchange field phase diagram of the bilayer is studied thoroughly. In Sec. IV, 
we first construct a Ginzburg-Landau functional to determine the transitions between the 
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different phases in the low interlayer coupling limit. Sec.V is devoted to the opposite limit 
of strong interlayer coupling. In conclusion, we discuss the conditions for the observation of 
field-induced superconductivity. 



II. ATOMIC THICKNESS SF/SF BILAYER 

We consider a superconducting ferromagnetic bilayer (see Fig.l) constituted of two super- 
conducting atomic layers, labeled as j — 1 and j = 2. In the whole article, we assume t <C E F 
where t is the interlayer coupling energy and E F the Fermi energy. As a consequence, Cooper 
pairs are localized within each plane. [15] Each layer j supports a superconducting singlet 
BCS coupling with the energy gap Aj and an internal exchange field hj. The Hamiltonian 
of the system can be written as 



H = E \ H i + H f cs + m / d2rA ?( r ) 



"t. (1) 

i=i,2 ' 

where A is the attractive BCS interaction constant and r is the two-dimensional coordinate 
within each layer. For the layer j the kinetic and Zeeman parts of the Hamiltonian are 
written together as 

^ = E^(p)*(p)^(p)> ( 2 ) 

p 

in which summation over repeated spin indexes a and a' is implied. Creation (resp. an- 
nihilation) operator of an electron with spin er and two-dimensional momentum p in the 
layer j is denoted Vv(p) ( res P- V^(p))- The exchange fields hj are assumed to be either 
equal (hi — h 2 — h) or opposite (h± — —h 2 — h). As a consequence the matrix ^ aa i is 
spin-diagonal, and the Zeeman effect manifests itself in breaking the spin degeneracy of the 
electronic energy levels according to 

W(p) = <W K(p) - ah j) > ( 3 ) 

where £(p) = p 2 /2m — E F . The s-wave singlet superconductivity is represented by the 
standard mean-field Hamiltonian 

h bcs = [A*(q)^ (p)^-t (-p) + h.c] , (4) 
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and the layers are coupled together by the hopping Hamiltonian 



H t = tJ2bPt(p)i>2a(l>)+h.C.} 



(5) 



p,cr 



In order to investigate the occurence of modulated superconducting phases (FFLO), we 
choose the following spatial dependence for the superconducting order parameter 



Ai(r) = Ae iq - r+i * /2 , A 2 (r) = Ae iq - r - 



(6) 



where q is the FFLO modulation wave vector and \ the superconducting phase difference 
between the layers. 

The above model can be solved exactly using the Green functions 



*5(p,p') = <^(p)i#t(p')> = sfr + jfiF+b), 

G jk (p, p') = - (ip^(p)^(p')) = 6(p - p' + q)G ifc (p), 



(7) 



where j and k are the layer's indexes. The brackets mean statistical averaging over grand- 
canonical distribution. [29] 

We obtain the following Gor'kov equations in the Fourier representation: 

/iw-£i T (p + q) -t A 1 \ /Gn(p + q)\ /A 

-i icu-^2t(p + q) A 2 
A* iw + 6i(p) t 



\ 



A* * ioj + U(p)J 



G 2 i(p + q) 
*£(p) 

v p 2i (p) y 








(8) 



where u = (2n + l)nT are the fermionic Matsubara frequencies. 

In quasi-2D superconductors [8, 30] the maximal FFLO modulation amplitude is of the 
order of (£o)~\ £o being the typical superconducting coherence length. This means that 
with a good approximation we can consider £jf(p + q) = £(p) — hj + vp.q, being the 
Fermi velocity vector in the plane. 

Solving the Gor'kov equations (8) yields the anomalous Gor'kov Green function for the 
j = 1 SF layer 



a 2 A* + t 2 A* 2 



where 



11 a x a 2 + (2A 2 cosx-/3)t 2 + t 4 ' 
aj = A 2 — ujj+ujj- and (3 = cui + lu2+ + uji-UJz-, 



(9) 



(10) 
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with Uj± = iu±^(p) — hj and ujj± = iuj ± £(p + q) — hj . Similar equation holds for F^. Note 
that in the case where a lattice made of such SF/SF bilayers is considered, the generalized 
anomalous Green function is obtained by replacing t 2 by \t + t'e~ %Pza \ 2 in Eq.(9), p z being 
the projection of the momentum p along the z axis and a the period of the lattice. As 
a consequence, a finite inter-bilayer coupling t' introduces an anisotropy in the dispersion 
relation which leads to a broadening of the electronic excitation levels. 

In the absence of tunneling t = we retrieve from Eq.(9) the anomalous Green function 
of a quasi-2D superconductor with the exchange field hi 

11 ai A* - (iu; + t n )(iu> - fry 1 ) 

Although the dependence on momentum has been removed for simplicity, notice that £^ = 
£jt(p + q) and ^| = Cj|(p). The set of basic equations (8) must be completed by the 
self-consistency equation 

/+oo 
diF±. (12) 
-°° 

Close to the critical temperature T c of the second-order phase transition, the order 
parameters Aj are small and Eq.(12) can be written as 

^. § /;>(«4 

where T c0 is the critical temperature for the 2D superconducting single layer in the absence 
of exchange field, namely for h — t — 0. At zero temperature, it is convenient to write 
Eq.(12) as 

where A = A(T = 0, h = 0, t = 0) is the superconducting order parameter at T = in the 
absence of exchange field and interlayer coupling. 



III. UNIFORM SUPERCONDUCTING STATES 



In this section, we investigate phases with uniform superconductivity within each layer. 
We obtain the critical temperature of the second-order superconducting (S) to normal metal 
(A) phase transition and the order parameter A(T, h, t) as a function of the temperature, the 
exchange field h and the interlayer coupling t. We consider both parallel (P) and antiparallel 



7 



(AP) orientations of the magnetizations, the superconducting phase difference being either 
X = or x — 71 ■ We also calculate the Josephson interlayer current when the bilayer is 
connected to external superconducting leads. Most of these results are obtained in the 
perturbative limit of small coupling between the layers t <C T c0 . The field is also assumed to 
be sufficiently small to prevent the occurence of a spatial modulation of the superconductivity 
within the planes. Study of nonuniform phases and strong coupling t ^> T c0 are respectively 
postponed to Sec. IV and Sec. V. 



A. Critical temperature 

We consider the second-order phase transition between the normal metal and the uniform 
BCS superconductor. Thus the order parameters A x and A 2 are small and the anomalous 
Gor'kov Green function (9) can be linearized in the following form 

F+ _ -^ 2+ £ 2 -A*+t 2 A* 

(Cl>i_CJ 2 - - t 2 )(uJi + UJ2 + ~ t 2 ) 

where £ J(T = £ — ahj. Similar equation may be found for F^ 2 . We first consider the 'parallel 
(P) orientation of the magnetizations, namely h — hi — hi. The first possibility is x = 
wherein the layers have the same superconducting order parameters A 1 = A 2 = A. In this 
situation the anomalous Green function obtained from Eq.(15) is denoted (F^) P '°. The 
identity 



Pfi 

c 



and Eq.(13) yield the following implicit equation for the critical temperature T c 

= - Re *(^w)' < 17 » 

where *f?(x) denotes the Euler digamma function. Therefore the interlayer coupling disap- 
pears from the self-consistency equation and Eq.(17) is identical to that for the 2D monolayer 
in a uniform exchange field: the bilayer is equivalent to a single layer in the neighborhood 
of superconducting to normal state transition. [2] The critical temperature of the bilayer 
decreases when the exchange field h increases. The equation (17) describes the second- 
order phase transition between the normal metal and the uniform superconductor which 
is realised only for fields smaller than the tricritical one h* = 1.07T c0 . For larger fields, 
superconductivity becomes nonuniform. 



A second possibility is the P orientation with x = tt phase difference between the layers. 
Now the anomalous Gor'kov Green function is denoted (F^j P,n . Then 



+00 



P,7T 7T 

2 



+ 



uj — i(h + t) lu — i(h — t) 



(18) 



and the self-consistency relation Eq.(13) yield a critical temperature T p ' w given by 



In 



c 



cO 



/ 1 \ 1 ( 1 h + at\ 

From this expression one may notice that superconductivity in the 7r state is destroyed by 
a combination of two effective exchange fields h ±t. In the small interlayer coupling limit 
t < T c , Eq.(19) becomes 

T/> , T p '° ( t V „ / 



ln ^= ln tr-feJ M^J- (20) 

where the function (x) is defined and represented in Appendix A. In the regime of low 
magnetic fields, namely for /i/27rT c p,7r < h* /2tcT* = 0.3, the factor K 3 (h/2%T ( P ' n ) is positive 
and thus the critical temperature is smaller in the 7r superconducting state than in the 
state. However the situation may be inverted if /i/27rT P ' ?r > 0.3. Moreover, along the 
critical line, the value h*/2irT* = 0.3 corresponds to the tricritical point, h* 1.07T c0 
and T* ss 0.56T c0 , where FFLO nonuniform states appear. As a consequence one expects 
competition between the 7r superconducting phase and FFLO phases in the neighborhood 
of the tricritical point. This competition will be detailed in Sec. IV. 

Let us focus on the case of AP orientation h — hi — —h 2 . Following the same procedure 
as previously, the equations for the critical temperatures T c Apx are obtained. In the limit 
t <T c0 , it reads 

rpAPfi _ rpPfi U2 

■ " 27rT c AP '°t 2 Y (21) 



± = -2nT^t 2 Y U — 2 (22) 

T po ^ 2 + ^ 2 ) 2 

for x = 7r. From Eqs. (21,22) the critical temperature is clearly higher in the phase than 
in the ir phase. Therefore the phase is the more stable in this region of the (T, h) phase 
diagram, i.e. in the vicinity of the critical temperature and for low fields h < h* = 1.07T c0 . 
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for x — 0, and 



In conclusion, the bilayer is always in the superconducting state for temperatures close 
to the critical temperature, whatever the relative orientation of magnetizations is. A spin- 
valve effect is also present: the critical temperature is higher for the AP orientation than 
for the P orientation of the magnetizations . 



B. Zero temperature superconducting gap 

For t — and low fields h < A /a/2, it is well-known that the zero temperature gap 
A(T = 0,/i, t — 0) = A is field- independent. [3] For small interlayer coupling t <C T c0 , the 
anomalous Green function (9) may be expanded to the second order in t as 

F+ 1 2 a x e ix - 2 A 2 cos x + P 
A\ ai a\a 2 

where the full nonlinear dependence on A is kept in a±,ai and [3. Then self-consistency 



relation (14) becomes 



4 = ^/fe^- 2A2 ) c ^^ (23 ) 
A J J 2tt a 2 a 2 v ; 



where A = A(T = 0,h,t) and x = or X = 71 ■ Using the preceding equation in the P 
orientation we obtain A p '° = A p,7r = A , either for or n phase difference. As a result, the 
superconducting gap A(T = 0,h,t) is not affected by a small interlayer coupling, at least 
at the order of t 2 . The superconducting condensation energy gain has also been calculated 
and the zero state found to be more stable than it state. 

For the AP orientation, the superconducting gap A^^ = A AP,0 (T = 0, h,t) is given by 



A AP,0 t 2 

In — - — = — 

A 2 



1 An - 2h 2 h 
+ , / a o , ono/o arcsm ■ 



A 2 ,-/* 2 ' /i(A 2 -/i 2 )3/2— A A _ ' (24) 
for zero phase difference. The unphysical divergence at h — > A is removed by terms of 
higher order in t. Expression (24) is the main result of this paragraph and reduces to 

, AAP '° 4t 2 h 2 

in the small field regime h A . Therefore in the x = state and for AP orientation, the 
order parameter is suppressed by the exchange field in the small coupling limit. This is sur- 
prising because AP orientation was expected to weaken the effective exchange field and thus 
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enhance superconducting properties. Nevertheless such a decrease of the superconducting 
order parameter has already been found in a ballistic atomic-scaled F/S/F trilayer.[22, 23] 

For the AP orientation and n phase difference, the gap A AP ' n (T = 0, h, t) = A is field 
and coupling independent. Moreover the energy of the x — 71 state does not depend on the 
relative orientation of the magnetizations. 

To summarize, the lowest energy corresponds to the (P, 0) phase. The (P, re) and (AP, n) 
phases are degenerate with a somewhat higher energy than the (P, 0) phase. Although 
we have not performed the energy calculation in the case where the magnetizations are 
antiparallel and the phase difference is 0, we believe that the highest energy corresponds to 
the (AP, 0) phase since its order parameter is the smallest one. 

C. Superconducting gap versus temperature: inversion of the proximity effect 

We now extend our investigation of the superconducting gap to finite temperatures. In 
order to determine the gap A(T,h,t) as a function of temperature T, exchange field h 
and coupling t, we analyse numerically the self-consistency relation (14) using the exact 
anomalous Gor'kov Green function (9). The result is shown schematically in Fig. 2. For 
the P orientation and x — 0, the superconducting gap A P,0 (T, h,t) is the same as the gap 
A(T, h, t — 0) of a single layer whereas for x = n the superconducting gap A p ' n (T, h,t) is 
lowered by finite interlayer coupling. For the AP orientation and x — , the gap is smaller 
than A(T,h,t = 0) for T < T; and larger for T > T { where the inversion temperature 
Tj = Ti(h) depends only on the exchange field in the small interlayer coupling limit (see 
Fig. 2 inset). This phenomenon has been called inversion of the proximity effect. [22, 23] 
Moreover, the gap A p '° is larger than A ,7r for all temperatures. 

According to these results, one may suggest several experiments. First we consider a 
bilayer with magnetizations pinned in the AP mutual orientation. By lowering the temper- 
ature, a 0-7T transition is expected at some temperature T n . In the small interlayer coupling 
limit, this temperature T w (h) is a function of the exchange field only (see Fig.2 inset). In 
contrast, the state is more favorable energetically for all temperatures in the case of mag- 
netizations pinned in the P orientation. As another illustration we consider samples where 
the relative orientation of magnetizations is free. Then the orientation is chosen by the 
system to minimize its energy. Cooling such a bilayer will result in a switching from the AP 
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orientation to the P orientation at the inversion temperature Ti(h). The same prediction 
was made recently in a ballistic F/S/F trilayer.[22, 23] 



D. Josephson current at T=0 

Here we consider that the SF/SF bilayer is connected to superconducting electrodes. In 
this set-up, one may impose an arbitrary superconducting phase difference \ between the 
SF layers, and thus a non dissipative Josephson current flows through the bilayer in the 
direction perpendicular to the planes. This interlayer Josephson current is evaluated here in 
the tunneling limit t <A and at zero temperature. Within the Green functions formalism, 
the general formula for the interlayer Josephson current is 

^XtKoJJ^^ (26) 

where N 2 d = m/(2nh 2 ) is the two-dimensional density of state per spin direction and unit 
surface. Solving exactly the Gor'kov equations (8) leads to 

t(^ 2+ -t 2 -A*A 2 ) 
21 aia 2 + (2A 2 cosx-/3)t 2 + t 4 ' 1 ' 

The function G\ 2 is obtained from Eq.(27) by permuting the layer indexes 1 < — > 2. The 

corresponding anharmonic current-phase relationship is given by 

2et 2 N 2D f[du 2A 2 sinx 



II' 



h J J 2ix a x a 2 + (2 A 2 cos x - I3)t 2 + t 4 

In the tunneling limit t « A , the interlayer Josephson current becomes sinusoidal as a 

function of the phase difference, 

• = • ff 2 A 2 sinx , . 

3 2vr ^(A 2 + e 2 + (^ + ^i))(A 2 + e 2 + (^ + ^/ i2 ))' 1 ' 

where jo — 2eN 2 Dt 2 /h . The second harmonic sin 2% has also been evaluated and is smaller 

than the first one by a factor (t/Ao) 2 . The preceding equation (28) yields the current-phase 

relation both for parallel (P) and antiparallel (AP) orientation of magnetizations. In the 

parallel case the critical current does not depend on the field as already reported in other 

systems since jp = jo sin [31, 32] For the antiparallel orientation and to the lowest order 

in t, the current-phase relation reads 

Jap = jo/i (£J sinx, (29) 
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where fi(x) = arcsinx/ (x\/l — x 2 ) . Therefore the critical current increases with the ex- 
change field h and even diverges for h = A = A AP (T = 0). Of course this divergence is 
unphysical and should disappear if all orders in t were taken into account. In Fig.3, the 
critical current is shown as a function of the exchange field both for P and AP orienta- 
tions. Recently, the issue of the Josephson coupling between two clean SF layers through an 
insulating layer was considered using Eilenberger equations [31] or Bogoliubov-de Gennes 
formalism. [32] Similar results as ours were obtained: the critical current increases with h 
only if three conditions are met: low temperature, very weak coupling between the SF lay- 
ers and AP orientation. Otherwise the presence of an exchange interaction suppresses the 
Josephson current. Using Usadel equations, Krivoruchko demonstrated that this statement 
holds in the diffusive regime for which the divergence for h = A is remplaced by a regular 
peak. [33] 

IV. PHASE DIAGRAM OF THE WEAKLY COUPLED SF/SF BILAYER 

From now on, we consider the SF/SF bilayer only for the parallel (P) orientation. Hence 
the results obtained in the next sections may be also applied to a superconducting bilayer 
in an external in-plane magnetic field. The present section is devoted to the weak coupling 
regime t « T c0 . In contrast to the low field restriction of Sec. Ill, regions of the phase 
diagram with h/{2^T) > 0.3 are also investigated here. Then competition between the 
FFLO and it phases is expected to take place. Particular attention is paid to the vicinity 
of the tricritical point given by h ~ 1.07T c0 and T « 0.56T c0 .[2] In order to examine this 
narrow region of the (T, h) plane, we construct a Ginzburg-Landau (GL) functional from 
the Gor'kov equations used in the previous sections. In the past Buzdin and Kachkachi [34] 
derived a generalized Ginzburg-Landau (GL) functional for a single SF layer that describes 
the FFLO superconducting state near the tricritical point. Here we extend this functional 
to a SF/SF bilayer for which it is possible to have not only FFLO modulation within the 
planes but also x = 71 superconducting phase difference between the planes. For x — 0, 
the physics of the bilayer is independent of the coupling and thus the Buzdin-Kachkachi 
GL functional is retrieved. In contrast for x = 7r, we obtain a free energy functional which 
depends on the interlayer coupling t and leads to the presence of a superconducting n phase. 
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A. Ginzburg-Landau free energies 



The free energy of the SF/SF bilayer in a uniform superconducting state with x — 
(U — state) is given by (see details in Appendix B) 



Fu- (A,f)=f 



+ b 



(30) 



with 



(31) 



where h = h/2irT and A = A/2nT are respectively the reduced exchange field and the 
reduced order parameter. Here we retrieve the well-known case of a single SF layer. For 
reduced exchange fields lower than the tricritical one h < h*, the transition between the 
superconducting and the normal states is a second-order one since e < 0. The critical line is 
given by the equation f(h, T) = 0. For higher fields, the transition becomes a first-order one 
because e > and b > 0. As for any first-order transition, two conditions must be fullfilled. 
On one hand, the free energy (30) is minimized, {dFu^/ dA) a=Ai = 0> anc ^ 011 the other 
hand the free energies of the superconducting and normal phases are equal, Fu-o(Ai) = 0. 
Hence in the (T, h) plane the equation for this first-order line is fi(h, T) = e 2 /(Ab) , and 



A l 



e/(26). It is 



the jump of the superconducting gap at the transition is given by 
well-known that this scenario is not realized because it is replaced by a transition between 
normal metal and nonuniform superconductivity. [8] Nevertheless, this first-order transition 
provides a useful energy scale A x and a reference line in the (T, h) plane that will be used to 
construct a universal phase diagram, namely a t-independent phase diagram valid for any 
weakly coupled SF/SF bilayers, see Sec. IV. C. 

The SF/SF bilayer may support opposite order parameters on the layers, superconduc- 
tivity being still uniform within each SF plane. In this so-called U — n state, the free energy 
of the bilayer depends on the reduced interlayer coupling t = t/2nT according to 



Fu_ n {A, f, t) = (r - 4eP + 8bt 4 ) 



(e - 12bt 2 ) 



A 



+ b 



A 



(32) 
(33) 
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For low reduced fields h < h*, F[/_ (A,r) < Fu^(K,f,t). Hence the uniform supercon- 
ducting phase with x = is more stable than the it phase, as already found in the Sec. III. 
Interestingly for higher reduced fields h > h*, this n phase is in competition with FFLO 
nonuniform superconducting phases having either \ = (FFLO — 0) or \ = 71 {FFLO — n). 

According to Buzdin and Kachkachi the order parameter A(x) = Acosgx leads to the 
lowest energy. [34] For the FFLO — phase, the corresponding free energy reads 



F LO -o(A, Q, f) = (^-2eQ 2 + 66Q 4 ) 



3 5ft 2 
8 e+ 16 Q 



+ — b 
16 



A 



(34) 
(35) 



whereas for the FFLO — n phase, the free energy depends on the interlayer coupling t in 
the following manner 



F LO -„(A, Q, f, t) = F LO - (A, Q, f) 

+ (-2et 2 + 4bt 4 + 24bt 2 Q 2 ) 
9,^ 



A 



(36) 
(37) 



The notation Q = VFq/ {4\f2nT) is introduced in Appendix B. 



B. Competition between FFLO-0 and U-7r phases 



Now we proceed to analyse the above free energies in order to determine the critical line 
between the normal and the superconducting states. We will also describe the nature of the 
various superconducting states and what kinds of transitions are encountered. We focus on 
the vicinity of the tricritical point. Then e = e(h — h*) is a linear function of the exchange 
field with e > 0, whereas b > is nearly field and temperature independent. 

For high reduced fields h > h*, it appears that the U — and the FFLO — n never lead 
to the highest critical temperature. Then we emphasize the competition between the two 
remaining phases, namely U — 7r and FFLO — 0. The minimization of the GL functional 
(34) leads to a modulation wave vector given by Q 2 = ej (6b) in the limit A — > 0, i.e. near 
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the critical line. The U — n phase is more stable than the FFLO — phase under the energy 
condition Fu-^A, f, t) < F LO -o(A, y/ej (66), r), or equivalently for 

2(3 - V3)bP < e < 2(3 + V3)bP. (38) 

Therefore the uniform ir phase is "inserted" within the usual FFLO superconducting 
state. The upper and lower values of h between which this new n phase is stable depend on 
the particular value of the coupling t. It is convenient to define a dimensionless generalized 
coordinate rj = e/(2bt 2 ) that quantifies the "distance" from the tricritical point along the 
S/N transition line. Indeed i] = at the tricritical point and the n phase is settled in the 
region (3 — \/3) < r\ < (3 + y/3). As shown on Fig. 4, going along the critical line from low 
to high fields, one expects the sequence of superconducting states: uniform in the planes 
with x — for rj < 0, FFLO modulation along the planes with x — for < rj < (3 — V3), 
then uniform n state for (3 — y/3) < rj < (3 + VS) and finally FFLO modulation along the 
planes with x — for i] > (3 + \/3). In all cases, the sign of the A 4 coefficient in the GL free 
energy is always positive and thus transitions between these superconducting states and the 
normal metal are second-order ones. 



C. Universal (r, 77) phase diagram 



We now construct the phase diagram around the tricritical point. Because each value of 
the coupling leads to different transition lines, we introduce the following mapping of the 
thermodynamic variables 

A f 

6 = j-,t=1 (39) 
Ai n 

in order to obtain a universal phase diagram valid in the small coupling regime. This 
mapping makes use of the energy scale Ai and of the function f\ related to the first-order 
transition between the normal state and the uniform superconducting state, see Sec. IV.A. 
Then the free energies for the uniform superconducting phases Eqs. (30,32) become 

Fu_ (5,t) 2 4 6 



and 



T|^-2ic+i<5r, ( 4 o) 



16 



Fu_„(6,T,ri) ( /l 1 



7 - S ( - " ^2 ) ) IJ ! 2 (ID 
6\ ,,,4 , i 



-2^1--j|*| 4 + |*r\ (42) 

where F = e 3 /(86 2 ). First it is straightforward to minimize Fu- (5, r) and Fu- n (S,r,T]) 
with respect to 5. Then replacing the reduced gap 5 by its equilibrium value, one obtains 
the equilibrium energies Fu_ q (t) and Fjj_ n (r,r]) of each superconducting phase. These 
energies are functions of both field and temperature via the dimensionless thermodynamical 
variables r and r). 

Using the same scaling Eq.(39) the free energies of the FFLO — phase, 



F LO . (5,Q,r) = fr_ _ 86 2 _ 246* ^ |JC|2 



and of the FFLO — n phase, 

FLo-n(6,Q,T,ri) F LO -o(5,Q,t) 



Fq F C 







/ 4 4 486^ 2 \ ... 2 

+ («) 

are also obtained for the order parameter A(x) = Acosgx. The equilibrium energies 
Flo-o( t ) and Flo—k(t,v) °f the modulated phases are obtained after minimization of 
Eqs. (43,45) with respect to <5 and Q .Then Eq.(43) enables to study the second-order phase 
transition between the normal metallic state and the nonuniform FFLO state. Under the 
assumption of second-order phase transition, it is sufficient to consider the GL free energy 
up to the S 2 order. Then the free energy (43) is minimal for Q 2 = ej (66). For this particular 
modulation, the critical FFLO/N line is given by r = 4/3. 

We now consider the transition lines between the various superconducting states obtained 
in the previous paragraph, in particular the FFLO — O/U — ir, the FFLO — 0/U — O and the 
U — ir/U — transitions. Let us focus on the transition between the uniform phases U — 
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and U — 71. Solving Fu-oij) = Fu- n (r,r]), we obtain the critical line t(t)) which corresponds 
to a first order U — 0/U — n phase transition. However, this transition is not realized (see 
Fig. 4) because the transition to the nonuniform superconducting state occurs before. 

The FFLO — O/U — transition line is obtained in a similar way. The equation Fu_ (r) = 
Flo-o{ t ) has the solution r « 0.913. At this value of r, the system undergoes a first order 
phase transition from the uniform state to the modulated FFLO state. Adding higher 
harmonics to the order parameter A(x) = Acosgx + A'cos3gx + gives a more accurate 
evaluation, namely r 0.859. 

Finally, the FFLO — 0/U — n transition line is obtained from F LO -o(t) = Fu^{t,ti) and 
shown in Fig. 4 in the (r, rj) plane. This transition is a first-order one. 

Using the mapping (T, h, t) — > (r, rj) of the thermodynamical variables, we have obtained 
a universal phase diagram Fig.4 of all weakly coupled SF/SF bilayers in the vicinity of the 
tricritical point. An important feature of this phase diagram is the presence of a supercon- 
ducting 7r-phase. As an example, the phase diagram has been redrawn in the (T, h) plane 
on Fig.5 for a particular value of the coupling t. 

V. PHASE DIAGRAM OF THE STRONGLY COUPLED SF/SF BILAYER 

Now we consider the ballistic SF/SF bilayer in the regime of strong interlayer coupling 
limit t ^> T c0 and low temperature. A very unusual 7r-superconducting state is found between 
a lower h\l2 —t — Aq/4£ and an upper hup —t + Aq/4£ critical exchange field, and below 
a maximal temperature of the order of T c 2 /t. Therefore field-induced superconductivity is 
obtained above hfj^ within the BCS theory of superconductivity. The underlying physical 
mechanism is the compensation of the Zeeman splitting by the energy splitting between 
bonding and antibonding electronic states of the bilayer, see Fig. 6. [27] Thus the new zero 
temperature paramagnetic limit hup = t + Aq/4£ may be tuned far above the usual one 
[8] h — A merely by increasing the interlayer coupling. This compensation also occurs 
for small coupling, but the 7r-superconducting state is then less energetically favorable than 
the usual O-superconducting phase as demonstrated in Sec. IV. Therefore the (T, h) phase 
diagrams are topologically distinct in the opposite limits of small (Sec. IV) and strong (Sec.V) 
coupling. We first analyse the second-order superconducting/normal phase transition in 
Sec.V.A. Then the first-order transition between uniform superconductivity and the normal 
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state is discussed in Sec.V.B. 



A. Second-order phase transition 



Here we study the second-order phase transition between the ir superconducting state and 
the normal metal state, as a function of the field. We start from the the linearized anomalous 
Green function (9) for arbitrary coupling t and n superconducting phase difference, 

A*, [t 2 -(iw + h- 2 ] • [t 2 ~(ioj + h + 2 ] ' 
From this equation and the self-consistency relation (12), the critical exchange field h is 
shown to satisfy 



h c + t + y/(h c + t)*-X* . h c -t+ y/(h c -t)*-X* = Ah 2 Q , (47) 

where X = \q\vp, and h = A /2 is the critical exchange field for the second-order super- 
conducting phase transition in a two-dimensional monolayer. One must then find the value 
of X which maximizes the critical field h c . If the n phase is assumed to be uniform inside 
each plane, namely if q = 0, Eq.(47) merely reduces to \h\ — 1 2 \ = h%. The lower and 
upper critical fields are respectively given by h c = t ± hl/2t, in the limit t 3> h . Thus 
at zero temperature and strong enough coupling, the superconductivity destruction follows 
a very special scenario. At low fields, superconductivity is first suppressed as usual at the 
paramagnetic limit /ifp = Ao leading to the normal metal phase. Then further increase 
of the field leads to a normal to superconducting phase transition at the lower critical field. 
This superconducting n phase is finally suppressed at the upper critical field. This is a new 
paramagnetic limit which may be tuned far above the usual one merely by choosing the 
coupling t greater than A . Thorough analysis of Eq.(47) shows that the upper critical field 
is even increased by an in-plane modulation in analogy with the case of the two-dimensional 
FFLO phase. [8] The upper critical field is maximal for the choice X = |q| vpjl = \h c — t\, 
and then Eq.(47) reduces to 



\h r . - 1\ . 



h c + t + 2\/h r t 



Ah 2 0: (48) 



that gives the upper and lower fields h^*] ow — t ± h^/t in the t > A limit. Note that the 
period of the modulated order parameter |q| _1 = £o(t/A ) is larger than the corresponding 
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period in the two-dimensional FFLO phase which coincides with the ballistic coherence 
length £ = v F /A . [8] 

Furthermore one may derive the full temperature-field phase diagram using Eqs. (13,46) 
and the result is shown in Fig. 7. When the temperature is increased, the lower critical field 
increases whereas the upper one decreases. Along the upper (resp. lower) critical line the 
FFLO modulation is lost at some temperature T* p (resp. Tj* ow ). For higher temperatures a 
uniform n phase (U — it) is recovered and the temperature dependence of the critical field 
is given by 



T c 2 ^ \2J \2 2txT J 

a=±l 



(49) 



where \1/ (x) is the Digamma function and ^(1/2) = — C — 2 In 2 ~ —1.963, C being the Euler 
constant. Finally the lower and upper critical lines merge at field h c = t and temperature 
Tm — 7re Td)/ (4£) in the limit t ^> T. Therefore the field-induced tt superconductivity is 
confined to temperatures lower than Tm- The structure of these U — it and the FFLO — tt 
phases is reminiscent of the corresponding U — and the FFLO — phases although the 
former are shifted to higher fields and lower temperatures than the later. 

Above results were obtained for relatively strong coupling. For lower coupling t ~ A , 
the U — it and the FFLO — n phases merge continuously into the usual x = phases as 
shown in Fig. 8, and finally disappear for t slightly smaller than A . From an experimental 
point of view, one might choose a system with intermediate coupling t small enough to settle 
the n phase island in an available range of temperatures but also large enough to separate 
the n phase island from the usual superconducting phases with x — 0. In the general SF 
multilayer case the inter-bilayer coupling constant t' needs to be sufficiently high to prevent 
from superconductivity destruction by 2D fluctuations but also sufficiently low to preserve 
the effect of field-induced superconductivity. [27, 35] 



B. First-order phase transition 

In the following we investigate the first-order U—n/N transition to determine whether it is 
more or less favorable than the above studied second-order transition. The zero-temperature 
superconducting order parameter A = A(t, h, T = 0) is calculated from the self-consistency 
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equation (14) for P orientation of magnetizations and x — 71 phase difference. At zero 
temperature, the difference between the energy Es of the superconducting state and the 
energy E N of the normal metal state is given by [29] 

A 



E S -E N = j 
o 



/ ' F+ /, , rx\ du) 



dd (50) 



In the limit t — > 0, we retrieve the well-known case of the single SF layer. [2, 3] Then 
the self-consistency relation (14) admits two branches of solutions. The lower branch A = 
a/A (2/i — A ), labelled (2) in the inset of Fig. 9, corresponds to a positive energy cost 
Es — En. Thus this superconducting solution is never realized. The actual superconducting 
gap is given by the upper horizontal branch A = A , (1) in the inset of Fig. 9, which 
corresponds to the energy difference 

E s -E N = ~(Al-2h 2 ). (51) 

Hence the superconducting phase is settled for low fields h < A /a/2 with a field-independent 
order parameter A = A . For higher fields h > A /\/2, the system is in the normal phase 
A = 0. Finally the zero temperature gap exhibits a jump at h — A /\^2 which reveals the 
first-order transition from the uniform superconducting phase to the normal phase. 

In the opposite limit of strong interlayer coupling, we have obtained in Sec.V.A. that 
field-induced superconductivity with x = 71 phase difference occurs for fields close to t and 
at low temperatures. From the self-consistency equation (14) one obtains several possible 
solutions for the zero-temperature superconducting gap A = A(t, h, T — 0) as a function of 
the exchange field h, see Fig.9. For relatively low fields h^ 10 = A < h < h_ and for high 
fields h > h + , the bilayer is in the normal phase A = 0. The limiting fields h± are solutions 
of 

h±-t\ 4 ft 2 -hl\ , x 

For intermediate fields ranging between and h + there are three superconducting branches. 
Two of them, (2') and (2") are never realized owing to their energy cost Es — E N > 0. The 
third branch (1') requires more detailed analysis. Namely, it is given by the equation 

aT 2 ^^ + 1 = ' (53) 
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and the corresponding energy cost is 



Es — E N — 



ttA 2 



2 



+ 



7l(h + t) 2 

2 



1 - 




2 



+ 



7l(h-t) 2 

2 



(54) 



Analysis of these equations reveals that E s — E N is negative for h\ ' w < h < h u J where 



h^l = Jt 2 - A2/(2 v / 2) and h!$ = Jt 2 + Ag/(2v / 2). Hence the SF bilayer undergoes first- 



order transition at h — h\ oul and h = h up . This scenario is quite similar than the one for 
t — 0, but with a smaller order parameter jump at the transition. Moreover there are two 



In order to generalize the above gap calculations to finite temperatures and determine 
the first-order S/N transition line, we have solved numerically together the self-consistency 
equation (13) and the condition Es — En = 0. The result is given in the inset of Fig. 7. 

Collecting results from Sec.V.A and B. we obtain the full (T, h) phase diagram for the 
field-induced n superconductivity. Note that this n superconductivity reproduces the struc- 
ture of the phase diagram in quasi-2D superconductors [8] although it is shifted to higher 
fields and lower temperatures. 

VI. CONCLUSION 

In this paper we have studied a periodic array of SF/SF bilayers in the limit of small cou- 
pling between the different bilayers. The corresponding Gor'kov equations have been solved 
exactly, taking into account both in-plane FFLO modulation and arbitrary superconducting 
phase difference between SF layers. The superconducting state with zero phase difference 
is always settled in the low field regime, h < T c0 for parallel (P) orientation of the magne- 
tizations. For antiparallel (AP) orientation, the n state predominates at low temperatures 
over the state which is settled in the neighborhood of the critical line. Consequently if the 
system is pinned in the antiparallel orientation, we predict a transition from the usual x — 
superconducting state to the n state by cooling. 

While the critical temperature is higher for the AP orientation, the zero temperature order 
parameter is larger for the P orientation. This results in a crossing temperature T^h) below 




first-order transitions, respectively at h\J w and hup instead of one at h — A /\/2 . 
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which the P orientation is more suitable for superconductivity. This temperature has been 
calculated as a function of the exchange field. In an experiment where the magnetizations 
might be easily reversed, one therefore expects a transition from the AP to the P orientation 
by cooling the system below this crossing temperature. 

In the low interlayer coupling limit, a Ginzburg-Landau functional has been derived from 
the exact expression of the anomalous Gor'kov Green function. As a main result, we have 
obtained a n superconducting state located in the vicinity of the tricritical point (h*,T*). 
Details of the bilayer phase diagram are obtained in this framework, including the first-order 
transition lines between superconducting phases. Since increasing the interlayer coupling 
enlarges the ir phase region, experimental observation of such details of the phase diagram 
requires the use of SF layers with large enough interlayer coupling, namely t ~ 0.1T c0 . 

Finally the case of even stronger interlayer coupling, namely t ^> T c0 has been also 
investigated. It appears that at low temperatures the n superconducting state is settled 
for exchange fields of the order of t, which are well above the Chandrasekhar-Clogston 
paramagnetic limit. Thus this new paramagnetic limit may be tuned by varying the interlayer 
coupling. In the present article we have reported the detailed structure of the phase diagram 
in this regime of high magnetic field. The first-order U — tt/N transition line is also derived. 

We expect that our results may be applicable to compounds like Bi 2 Sr 2 CaCu208 un- 
der a magnetic field. Indeed such perovskite superconductors comprise tightly coupled 
superconducting CuO planes separated by BiO layers. However observing the field-induced 
superconductivity in a reasonable range of magnetic field requires relatively low critical 
temperatures which are realized in the heavily doped or underdoped regimes. Finally the 
latter effect is solely related to the compensation of the energy shift in the two layers sys- 
tems by the Zeeman splitting. So it should be quite general and might appear also in two 
band superconductors or in weakly coupled superconducting grains. Note that the inhomo- 
geneous superconductivity has been obtained in the absence of magnetic field in two-bands 
superconductors. [28] However since the n state was not considered in this latter work no 
field-induced superconductivity had been noticed. 

Bulaevskii [8] studied thoroughly Josephson coupling in periodic layered structures with 
one SF plane as unit cell. Here we have demonstrated that systems with several SF planes 
as unit cell exhibit qualitatively new phenomena like field-induced superconductivity. The 
simplest case, two planes per unit cell, has been studied here. It may be regarded as a 
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basic approach to understand the properties of more complex ferromagnetic superconducting 
compounds or artificial heterojunctions. 

We thank M. Daumens, M. Faure, M. Houzet and M. Kulic for useful discussions and 
comments. This work was supported, in part, by ESF "Pi-shift" Program. 

VII. APPENDIX 

A. Definition and properties of the functions K^(h) 



We define the function K\{h) by 



K 1 (h) = J2 Re 



n=0 



1 



1 



n+\ + ih n + \ 



(55) 



and for any integrer fi> 2 the function K^{h) is given by 

1 



o (n + | + ihj 



(56) 



The variations of K 1 , K 3 , K 5 , with h are represented in Fig.10. One can notice that in 
the vicinity of the tricritical point, i.e. h h* ~ 0.3, the functions Ki{h) and K 5 (h) are 
negative and of the order of unity. K 3 (h) cancels exactly at h = h* and becomes negative 
in the domain h > h*, which is studied Sec. IV. 



B. Ginzburg-Landau functional 

This part of the Appendix refers to Sec. IV of the paper. In the Ginzburg-Landau theory, 
the free energy is expanded in terms of the gap A, i.e. the order parameter, assuming the 
temperature close to T c . Originaly it was introduced as a phenomenological theory for 
superconductivity before the BCS theory. Here we derive the Ginzburg-Landau free energy 
from the full microscopic knowledge of our model in order to analyze the vicinity of the 
tricritical point. To do this, we consider the simplest case where the FFLO gap modulation 
is exponential, namely A(x) = Ae zqx , q being the in-plane modulation wave vector. It is 
known that this modulation structure is not realized to the benefit of the cosine modulation 
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discussed in the article's body. However, In Sec. II, the Gor'kov Green functions of the 
SF/SF bilayer were derived for a modulated order parameter A(x) = Ae' iqx and x = or 7r. 
This modulation structure is then convenient to calculate the coefficients of the generalized 
GL functional because the exact expression of the anomalous Green function (see Eq.(9)) 
is valid for this gap modulation structure, whereas it is unknown with the cosine structure. 
We first expand the exact anomalous Green function (9) and the self-consistency relation 
in powers of the gap A and the FFLO wave vector q. Then this self-consistency relation is 
interpreted as the stationarity condition for the Ginzburg-Landau free energy, which allows 
(by identification) to determine the coefficient of every term of the GL functional. 

In the x — case, the expansion of the anomalous Green function reads 

1 °° (— l) n+1 |A| 2n ~ 2 

# = 2^-a(t,q)r^ + a(-t,0)) n + ( ' ~ _t) ' (5?) 

where a(t, q) = iuj — h + t — vp.q. We first consider the case of uniform superconductivity, 
i.e. q = 0. After integration over £, we obtain: 

+oo 



/ 



, A i2n 



§-§ = s 9 n(u)J2(-lT b f ( , 



( 58 ) 



with b n = Jj222n = r ol / 2 )n \ ■ Note that the interlayer coupling t has disappeared in Eq.(58). 
We are now able to write down the self-consistency equation (13) as an expansion in powers 
of A 

T ~ _ _ 2n 

In— S^(-l) n b n K 2n +i{h) A =0, (59) 

where the functions are those defined in Appendix A, and A = A/2irT. This self- 
consistency relation may be interpreted as the stationnary condition 

^ = (60) 
<9A 

for the Ginzburg-Landau free energy with uniform order parameter within each supercon- 
ducting plane and x = phase difference between the planes. Close to the tricritical point, 
A is small and it is enough to retain only the first term in this infinite expansion as 

l n ^- boKl (h) +bl K 3 (h) A 2 -b 2 K 5 (h) A 4 = 0, (61) 



cO 
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where h = h/2ixT. By identification with Eq.(60) we obtain the GL free energy for the 
U — phase as a function of the variational parameter A and the thermodynamical variable 
h = V(2vrT): 



\n—-b Q K l {h) 

J-cO 



In— -K x {h) 

J-cO 



A 

+ b 1 K 3 (h)—-b 2 K 5 (h) 



A 





A 


2 ^) 


A 


4 ^ 5 (/i) 


A 






4 




8 





(62) 



which corresponds to Eq.(30). Note that in this usual state the same coefficients have 
been already reported in Ref.[34] 

The same procedure may be followed when the phase difference is ir. The anomalous 
Green function is then 

(-1)' 



# = 2^(e-a(U)r(e + aftO)) 



(63) 



and leads to the self-consistency relation which contains explicitely the coupling t, via the 
normalized coupling i = t/2nT: 



T 



c0 n=0 



K 2n+1 (h + t)+K 2n+1 (h-t) 



2n 



0. 



(64) 



From the latter expression we deduce that the coefficients of the GL free energy for the 
X = 7r state can be directly obtained using the coefficient of F[/_ (A) in which we replace 
K2n+i{h) by (K 2n +i(h + t ) + K 2n +i{h — t))/2. Finally the free energy of the U — n state is 



In 



T K^h + t) + K^h - t) 



+ 



T c o 2 
K 3 (h + t)+K 3 (h-t) K * 

A 



K 5 (h + t)+K 5 (h-t) 
16 



(65) 



which yields Eq.(32) in the small interlayer coupling limit t h. 

We have developped in a similar way the GL free energy in the case where the order 
parameter is modulated within each superconducting plane. Using the expressions (57) and 
(63) for the anomalous Green function of the bilayer with A(x) = Ae iqx FFLO modulation 
respectively in the x — an d X — 77 cases, one obtains the expansion of the self-consistency 
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equation in powers of A and of the FFLO wave vector q. Finally, after averaging over all 
possible orientations of the FFLO modulation vector, the self-consistency equation reads: 



In 



T 



- cO 



EE*- 1 ' 

n=0p=0 



n+p 



c n,pK2(n+p)+l{h) 



(VFq\ 
\4nTJ 



2p 







(66) 



for x = 0. The coefficients c„ iP are symmetric with respect to the expansion indexes n and p 

T(n + p+ 1/2) (n + p)\ 



Cn,p 



r(l/2) (p\) 2 (n\f 
and related to the coefficients b n by b n = c„ i0 . 

From Eq.(66) the GL free energy can be constructed using the method described in the 
previous paragraph for uniform phases. We retrieve all the coefficients already obtained by 
Buzdin and Kachkachi,[34] including the coefficients of the gradient terms of the generalized 
functional. To derive the free energy of the FFLO — phase, we have therefore used the BK 
functional with the cosine modulation which is effectively realized in each superconducting 
layer. As a result, it reads 

2 



F LO _ (A, Q, t) = g - 2eQ 2 + 66Q 4 ) 



8 16 V 



+ — b 
16 



(68) 



where Q = VFq/(Ay / 2TiT). In the x — 71 state the free energy has been derived from the BK 
functional in which the replacement 

K2(n+p)+l(h + t) + K 2 ( n +p)+l(h - t) 



K- 



2(n+p)+l 



th) 



(69) 



has been done in order to obtain the modified coefficients. Finally the free energy of the 
FFLO — 7r phase can be written as 



F LO _ ff (A,Q,r) 



f - AeP + 8bt 4 



2 (e - 126?) Q 2 + QbQ' 



+ T6 6 



(70) 



In the article body more convenient forms of Eqs. (68,70) involving the reduced quantities r, 
5 and r\ are used in order to derive the universal phase diagram. 
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VIII. FIGURE CAPTIONS 

FIG.l: SF/SF bilayer. The interlayer coupling constant is denoted t. The exchange 
fields hj can be either equal (parallel orientation) or opposite (antiparallel orientation). 
The superconducting phase difference between Ai and A 2 can be either (Ai = A 2 ) or it 
(A: = -A 2 ). 

FIG. 2: Schematic representation of the superconducting gap as a function of temperature 
for P orientation and x — (thicker solid line), AP orientation and x — (intermediate 
thickness line), and AP orientation and x = 71 (thiner line). All curves are given for the 
same value of the exchange field that is smaller T c (in energy units). The inversion of the 
proximity effect occurs at the temperature Ti(h) and the transition from state to ir state in 
the AP orientation at the temperature T n (h). Temperatures T^h) and T n (h) as a function 
of the field are shown in the inset. 

FIG. 3: Enhancement of the critical current j c with the field in the AP orientation (dashed 
line). In the P orientation j c does not depend on the exchange field h (solid line). 

FIG. 4: Universal phase diagram for weakly coupled bilayers, in (r, 77) coordinates. The 
critical line (solid line) corresponds either to a U — ir/N or to a FFLO — 0/N transition 
depending on 77. The transition between the nonuniform FFLO — superconducting phase 
and uniform U — (resp. U — n) phase is represented with dash-dotted (resp. dashed) line. 
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FIG. 5: Phase diagram in (T, h) coordinates, for t/ (2%T c0 ) = 0.07. Only the neighborhood 
of the tricritical point is represented. The lines have the same meaning than in Fig.4. Note 
that the n phase is settled in a very narrow region of the phase diagram. 

FIG. 6: Excitation spectrum. Usual singlet pairing (thin line circles) between opposite- 
spin electrons occupying the same orbital is affected by Zeeman effect. In contrast, is coupling 
(thick line) between two electrons occupying a bonding and an antibonding orbitals may lead 
to the cancellation of the Zeeman splitting. 

FIG. 7: Phase diagram for t = 3A ps 5.3T c0 . Thick (resp. thin) solid lines represents 
second-order transition between U — x (resp. FFLO — x) and normal metal phase (N) for 
X = and it. We expect the U — xi 'FFLO — x transition lines (not calculated) to be in the 
vicinity of the (virtual) first order U — x/-/V lines (dash-dotted). 

FIG. 8: Phase diagram for t = A ~ 1.76T c0 . All lines have the same meaning than in 
Fig.7. 

FIG. 9: Order parameter at T = in the (P, ir) state for t = 3A as a function of the 
exchange field (thick solid line). 

FIG. 10: Functions Ki(h), K 3 (h) and K 5 (h) defined in Appendix B. 
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